Late-time tails of a Yang-Mills field 
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We study the late-time behavior of spherically symmetric solutions of the Yang-Mills equations 
on Minkowski and Schwarzschild backgrounds. Using nonlinear perturbation theory we show in 
both cases that solutions having smooth compactly supported initial data posses tails which decay 
as t~ 4 at timelike infinity. Moreover, for small initial data on Minkowski background we derive the 
third-order formula for the amplitude of the tail and confirm numerically its accuracy. 
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I. INTRODUCTION 

^ n 

In a classical paper [1| Eardley and Moncrief proved that solutions of the Yang-Mills equations on the 3+1 Minkowski 

spacetime starting from smooth initial data remain smooth for all future times. A different proof allowing for initial 
data with only finite energy was given later by Klainerman and Machedon . Once global regularity was established, 
the problem of asymptotic behavior of solutions for t — > oo was studied by many authors [1, |j, H, [y, @] who obtained 
!— i ' various decay estimates using different techniques and assumptions about initial data. In this paper we are concerned 
with the simplest possible situation, namely spherically symmetric initial data with compact support. In this case it 
^ ' follows from the conformal method of Christodoulou that the Yang-Mills curvature decays as t~ 4 at timelike infinity 
5— i , The purpose of this paper is threefold. First, we rederive Christodoulou's result using the nonlinear perturbation 

, theory. The advantage of our approach lies in its wide applicability; in contrast to the conformal method which is 

| very powerful (in the sense of giving sharp decay rates) only for conformally invariant equations. 
CNJ ■ Second, we go beyond qualitative decay estimates and give the third-order formula for the amplitude of solution 
which provides a precise quantitative information about the tail. We wish to point out that although our result depends 
crucially on spherical symmetry, the assumption of compact support for initial data is made for simplicity and can 
be relaxed by imposing a suitable fall-off condition at spatial infinity (which can be implemented via appropriately 
weighted norms). However, some kind of localization condition is necessary in order to avoid a situation where the 
tail in time is induced entirely by the tail of initial data at spatial infinity (due, for instance, to nonzero charge). 

Third, we argue that the same tail is present in the scattering of spherically symmetric Yang-Mills fields off the 
Schwarzschild black hole. In this case the global existence of solutions follows from the work of Chrusciel and Shatah 
Q who generalized the proof of Eardley and Moncrief to arbitrary globally hyperbolic Lorentzian 4-manifolds. The 
late-time tail of the Yang-Mills field on the Schwarzschild background was studied in 0, however the fall-off t~ 5 
derived there on the basis of the linear perturbation analysis is not correct. As we shall see, the error in || is due 
to the fact that the late-time tail is not governed by the linearized evolution. At first sight that might seem odd 
but upon reflection it is easy to understand. A rough intuitive explanation is that the tail is a far-field effect hence 
the flat space tail t~ A is expected to persist in any asymptotically flat spacetime as long as the backscattering on 
the curvature does not produce a more slowly decaying tail. A similar example of the failure of linear perturbation 
analysis was recently observed in the scattering of skyrmions (To| . 
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II. MINKOWSKI BACKGROUND. 



We consider the Yang-Mills theory with the gauge group SU (2) and assume the spherically symmetric ansatz for 
the connection [ll| 

A = w T\dB + (cot 9 T3 + w T2) sin 9 dip , (1) 

where w = w(t, r) and r, (i — 1, 2, 3) are the usual generators of su(2). The Yang-Mills equations d * F = 0, where 
F = dA + A A A is the Yang-Mills curvature, reduce then to the semilinear radial wave equation 

to - w" - -jw(l - w 2 ) = 0, (2) 
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where primes and dots denote derivatives with respect to r and t, respectively. For our purposes it is convenient to 
define the function f(t, r) = (w(t, r) — l)/r and rewrite equation in the following form 

CI := / - /" - \l' + ^f = -f - \f . (3) 

Note that C is the radial wave operator for the I = 1 spherical harmonic. 

We consider the late-time evolution of solutions of equation ([3|) for smooth compactly supported initial data 

/(0,r) = ea(r), /(0,r) = e/?(r) . (4) 

The prefactor e is added for convenience and to emphasize that our initial data are assumed to be small. Regularity 
at the origin is ensured by the boundary condition f(t, r) ~ b(t)r for r — > 0. As follows from [3j such solutions decay 
to zero on any compact region of space as t — ► oo. To determine the asymptotic behavior of solutions we define the 
perturbative expansion 

/-e/i+£ 2 /2+e 3 / 3 + ..., (5) 

where ef± satisfies initial data (J3J) and all /„ with n > 1 have zero initial data. Substituting the expansion |[5J into 
equation ([3]) up to the third order we get 

Ch = 0, (6) 

Ch = --A 2 , (7) 
r 

Cfa = -A 3 --/i/ 2 . (8) 

r 

We solve these equations recursively. The first order solution is given by the general regular solution of the free radial 
wave equation for the I = 1 spherical harmonic 

ai(t-r) + a'{t + r) a(t - r) - a(t + r) 
fi(t,r) = + - 2 , (9) 

where the function a(£) is determined by the initial data 



DC.- 



«(0 =-\tJ a ( s ) ds + f)^)ds . (10) 

€ C 

For compactly supported initial data the function a(£) has compact support as well (note that the functions a(s) and 
P(s) in (fT0|) are odd extensions of initial data to the whole line), hence /i has no tail in agreement with Huygens' 
principle. 

To solve equations for the higher order perturbations we use the retarded Green's function of the operator C 

2 , „/2 _ /. _ .A2 

G(t - t', r,r') = [\r-r'\<t-t' <r + r'} ^ '- . (11) 

It follows from flT]) that the solution of the inhomogeneous equation Cf = N(t, r) with zero initial data has the form 
(using null coordinates u = t' — r' ,v = t' + r') [3(] 

t+r t—r 

/(*' r ) = ^2 j dv J K{t,r- : u,v)N(u,v)du, (12) 

\t—r\ — v 

where the kernel K(t, r; u, v) — (v — t)(t — u) + r 2 . In the second order, i.e. for equation l[Tj). the representation (jX3J) 
yields 

t+r t—r 

/ 2 (t,r) = -^ J dv f K(t,r;u,v)^^du. (13) 

\t—r\ — « 
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Somewhat surprisingly, Huygens' property is preserved in the second order. To see this, let us assume that a(£) = 
f° r |£| > R- Then, for t > r + R we may change the order of integration in (fT3"|) and rewrite it as (see Fig. 1) 



-i? t-r 



Performing the inner integral we get 



/ 2 (*,r) = 8r 



a(w) e? / a(it) 



(14) 



-ii 



(t — u) 2 — r 2 du \ (t — u) 2 — r 



du . 



(15) 



which after integration gives zero. Thus, f2{t,r) vanishes identically for t > r + R and consequently there is no tail 
up to the second order. 




FIG. 1: An illustration of the situation in equations (|14[) , (|16l) . and (|17[) . The observation point is located at (t, r) where 
t > r + R. The integration range is given by the intersection of the two shaded regions which depict the domain of dependence 
of the observation point and the support of the solution fi(t',r'). 



In the third order, i.e. for equation ([5]), the representation p^|) gives f% = /g + where 



fi 1] (t,r) 



1 

8^2 



t+r t—r 



dv / K(t,r;u,v)f 1 (u,v)du, 



2r 



\t—r\ -v 

\ Tdv7K(t 7 r;u,v) h ^ v)h ^ v) du. 



v — u 



(16) 
(17) 



t—r —V 



To calculate f^\t, r) for t > r + R, as above we change the order of integration and perform the integral over v with 
the result (using the abbreviation z = (t — u) 2 — r 2 ) 

R 



/ 3 (1) (t,r)=4r 



a{u)a'{u) 2 ( 4(t-u)a'(u)a 2 {u) 4((t - u) 2 + \r 2 )a z (u) 



+ 



du . 



(18) 



which has the following asymptotic behavior near timelike infinity (r = const and t — > oo) 



(t, r) ~ c\ r t 4 , ci = 4 / a(u)a \u) A du 



(19) 



In the formula above we replaced i? by oo in the limits of integration to emphasize that the result holds not only for 
strictly compactly supported initial data but also for initial data which fall off sufficiently fast at spatial infinity. 
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(2) 

To calculate the contribution to the tail coming from /g ' (t, r) we need to know both the leading and the subleading 
terms in the asymptotic expansion of f2{u, v) near null infinity (u = const and w — > oo). This calculation is deferred 
to the appendix where we show that near null infinity 

ft'(tQ , 2ft(u) 2g(u) 

h{u,v) = + ? + - ), (20) 

v — u (v — u) z (v — u) z 

where h(u) and g(u) are defined by (|32|) and l|33|) . respectively. Note that the first two terms in (f20f represent the 

"free" part of the iterate J2\t,r)\ as we shall see in a moment this part does not affect the behavior of / 3 (t, r) at 
timclike infinity. Substituting (|20p into (TT7|) and proceeding along the same lines as in the derivation of the expression 
(1191) we obtain the following asymptotic behavior near timelike infinity 



/f (t,r)~c 2 rt-\ 



C2 



— (h(u)a(u)) + g(u)a'(u) 
du 



du = -12 / a(u)a'(u) 2 du, (21 



where the last expression follows from (j3"3")l and integration by parts. Putting equations (fT!)]) and (j2~lj) together we 
finally get the leading asymptotic behavior near timclike infinity 



+oo 

\2, 



f 3 (t,r) ~ crT 4 , c = -8 / a(u)a'(u) 2 du . (22) 



This is our main result. We claim that the expression (|22|) provides a very good approximation of the tail for solutions 
having sufficiently small initial data. More precisely, we conjecture that for any given smooth compactly supported 
functions a(r) and (3{r) one can choose e such that for each fixed r > and t — > oo the remainder \t 4 f(t,r) — e 3 cr\ 
is as small as one pleases. The numerical evidence supporting this conjecture is shown in Fig. 2. The obvious issue 
remains as to whether the perturbation expansion corresponding to given initial data is convergent for sufficiently 
small values of e. Without a proof of convergence our analysis is not mathematically rigorous, nevertheless we believe 
that it gives a rather convincing and, most importantly, quantitative description of the late-time tail. 

A few remarks are in order: 

Remark 1: It should be clear from the above derivation that the simplicity of the final result (j22f is due to some 
amazing cancelations (notably those occurring in equations (|15[) and (|21[) 1 which in turn are attributed to the particular 
form of the nonlinearity of the Yang-Mills equations. In this respect the Yang-Mills equations are exceptional and 
particularly interesting mathematically; for most other nonlinear perturbations of the wave equation the tail is a 
second-order phenomenon which is much easier to analyze (e.g., see (To}). 

Remark 2: Note that all iterates fk behave as 0[(v — w) 1 ] near null infinity and therefore at each order of the 
perturbation expansion the sources behave as 0[(v — u) ]. For such sources one might expect by dimensional 
arguments that there would be a t~ 3 tail. Fortunately, due to the identity 



t+r 

(v - t)(t - u) + r 2 



(v — u) 3 



dv = (23) 



all coefficients of hypothetical t~ 3 tails vanish identically and thus all higher-order terms in the perturbation expansion 
decay as i -4 . This fact is crucial since otherwise the third-order approximation would break down for late times; for 
instance a nonzero fourth-order term ~ e 4 /t 3 would make the formula (|22|) useless for times t>l/e. 

Remark 3: Note that equation © has the scaling symmetry: if f(t,r) is a solution, so is f\(t,r) := Xf(Xt,Xr). 
Under this scaling the energy scales as E(f\) — XE(f), hence given any finite energy initial datum one can scale it 
down to an arbitrarily small amplitude and energy. Note, however, that for compactly supported initial data such a 
rescaling spreads the support by a factor 1/A and for this reason it cannot make large data smaller in the sense of 
our perturbation expansion. This follows immediately from the fact that all iterates fk(t,r) in {5} scale in the same 
way. In other words, the rescaling does not change the convergence properties of the perturbation expansion. 
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We plot (on log- log scale) the numerical solution fit, 

_2\ at \ In 2n _ / 2\ 



of the initial value problem (O-© for air) - 

.2 r i _3_ 



FIG. 2: 

rexp(— r 2 ), /3(r) = r(2 — r' 1 ) exp(— r A ) and e = 0.1, and compare it with the third-order approximation efi + + e a /3 
(produced by solving numerically the perturbation equations (6)- (8)). These two functions are indistinguishable at the scale 
of the figure so both are depicted by the single solid line. The contributions of the individual iterates are superimposed 
to demonstrate that the tail comes from f$. The fourth-order iterate fi serves as an estimation of the error. Note that 
fi has the same late-time slope (decay rate) as /a, in agreement with Remark 2. The fit of the function Ct -1 exp(A/t) to 
the full solution f(t,l) for late times gives 7 ~ 4.001 and C ~ 2.391 • 10~ which is ca. 4% off the third-order prediction 
e 3 c = e 3 155^37r/20736 w 2.295-10" 5 obtained by evaluating the integral in (|22]l. 



III. SCHWARZSCHILD BACKGROUND 



On the exterior Schwarzschild spacctime 



lls 2 f 1 dr" 

r 



1 - 



2m 



dr 2 +r 2 (d6 2 + 



sin 



r > 2m , 



the spherically symmetric Yang-Mills equation corresponding to the ansatz (JT|) takes the form 

-1 



1 



2m 
r 



1 - 



2m 



- \w{\ - w 2 ) = 0. 



When m = this equation reduces of course to ([2]) . In terms of the new variables 

x = r + 2m In ( 1 ] , hit, x) = wit, rix)) — 1 

V 2m / 

equation |25|) becomes 



•• d 2 h ( 2m\ 2 , 
h- —5 + I 1 ] — h 

dx 



!-^U(3/ i 2 + / l 3 ) 



(24) 



(25) 



(26) 



(27) 



where r — r(x). Dropping the nonlinear terms on the right side of (|27|l one gets the linear 1 + 1 dimensional wave 
equation on the real axis —00 < x < 00 with the effective potential V(x) — 2/r 2 — Am/r 3 . This equation describes the 
propagation of the dipole (I = 1) electromagnetic perturbation of the Schwarzschild black hole. For intermediate times 
the linearized approximation is very good; this stage of evolution has the form of exponentially damped oscillations 
dominated by the fundamental (i.e., least damped) quasinormal mode. We recall that quasinormal modes are solutions 
of the linearized equation satisfying the outgoing wave conditions h(t, x) ~ e^M*^) f or x _> ±00. In the case at 
hand the fundamental quasinormal mode has the eigenvalue k = 0.49653 — 0.18498? (in units where 2m = 1) [Hj]. 
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The quasinormal mode decays exponentially so for late times it becomes negligible and eventually a polynomial 
tail is uncovered. Since the pioneering work of Price [131 ] it has been known that the tail of the l-th multipole decays 
as t~ 2l ~ 3 , thus for the dipole the linearized theory predicts the tail t~ 5 and this is exactly the result derived in [fj]. 
We claim that this prediction is incorrect and the actual tail behaves in the same manner as in Minkowski spacetime, 
namely it decays as t~ 4 . Regarding equation ([25]) as the perturbation (for r ^> 2m) of equation ((2|), one can see from 
dimensional considerations that the failure of linearization is due to the fact that the linear terms in ([25]) corresponding 
to nonzero curvature (proportional to m) are of shorter range (using PDE jargon) than the nonlinear terms. Thus, 
the presence of the black hole should not alter the flat space tail t~ 4 . The numerical substantiation of this handwaving 
argument is shown in Fig. 3. 
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FIG. 3: Scattering of the Yang-Mills wave off the Schwarzschild black hole (with 2m = 1). We plot (on log-log scale) 
the numerical solution h(t,x = 1.5) of equation (|27l) for the initial data of the form of the "ingoing" gaussian h(0, x) = 
Aexp( — (x — xo) 2 1 s 2 ) with A — 0.85, xo — 3, s = 1.5. Fitting the exponentially damped sinusoid Q(t) = Be~ Tt sin(f2t + 8) on 
the interval 30 < t < 60 we get V = 0.184 and = 0.495, in perfect agreement with the known parameters of the fundamental 
quasinormal mode. The fit of the power law decay P(t) = Ct ' exp(D/t + E/t 2 ) for times t > 300 gives 7 = 3.9997. The sum 
\Q(t)+P(t)\ (depicted by the dashed line) provides a remarkably good approximation of the full solution for all t > 20. It should 
be pointed out, however, that our initial data were tuned a bit to maximize the effect of the nonlinearity. If the subdominant 
t~ 5 tail coming from the potential has a large coefficient, i.e. the tail behaves as Ct~ 4 + Ct~ 5 with C«C', then one has to 
wait for a long time before the true asymptotic behavior sets in (which might be misleading without an analytic insight). 



Unfortunately, for the Schwarzschild background we were not able to derive a quantitative formula, analogous to 
(|22[) . relating the amplitude of the tail to initial data. An attempt to repeat the perturbation analysis from section 11 
encounters serious difficulties on Schwarzschild background which are caused by the violation of Huygens' principle 
in 1 + 1 dimensions and the presence of the potential. It would be interes ting to pursue this problem further, perhaps 
borrowing ideas from an approach proposed some time ago by Barack Although Barack considered only the 

linear wave equation, we wish to emphasize that there are many similarities between his work and our analysis in 
section II. 



IV. CONCLUSIONS 



Using third order nonlinear perturbation theory we determined the late-time tail of spherically symmetric Yang-Mills 
equations on Minkowski background. We also gave heuristic arguments that the same tail is present on Schwarzschild 
background. In both cases we provided numerical evidence supporting our results. We hope that our approach will 
trigger more rigorous mathematical analyses of these physically important phenomena. 

We remark that the ideas presented here can be applied to other nonlinear wave equations. For example, one can 
show by similar methods that for the semilinear wave equation g^V M V ' v <p + \<fi\ p = on Minkowski background the 
tail decays as t x ~ v for p > 1 + v2j while on the Schwarzschild background the tail changes its character at p — 4 from 
linear (Price's law ~ t~ 3 for p > 4 15]) to nonlinear (<j> ~ t l ~ p for 1 + \/2 < p < 4). A systematic analysis of the 
competition between linear and nonlinear effects in scattering for semilinear wave equations in Minkowski spacetime 
will be given elsewhere [lfjj . 
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Appendix 



We derive here the equation ([20)) . Our starting point is the equation (|T3|) in which we relabel coordinates (it, v) — > 
(u', v') and use the retarded time u = t — r: 

u+2r u 

3 f , f u' + v' -2u+(uu' + uv' -u'v' -u 2 )/r „ , , , 

f 2 (u,r)=-— / dv' / ^ '-l-ft(u',v')du'. (28 

Ar J J v — u 

\u\ —v' 

We let e = 1/r and define the quantity 

u+2/e u 

3 f , , f u' + v' - 2u + e(uu' + uv' - u'v' - u 2 ) 2 , , , 
J(«,e):=-- / rfw'y ^— > -fl{u' 1 v')du' . (29) 

|m| — «' 

Expanding this in Taylor's series J(«, e) = A(u) + B(u)e + 0(e 2 ) we obtain 

oo u 

3 f , , f u' + v' -2u 



A(u) := l(u,0) = ~ I dv' I " y_ u "" tf(u',v')du', (30) 



|u| — u' 



5J(u,e) 



:= _:= fl (u) + h(u) , (31) 



£ = 



where 



7 / \ 3 f , , f uu' + uv' - u'v' - u 2 2 , , , 

= -g / du / 7/3^ fl{u\v')du' (32) 

|u| — 

c?(u) = 3^ a'(u') 2 du' . (33) 

— oo 

An elementary calculus exercise shows that 

h'(u)=A(u). (34) 
Putting all the above equations together and noting that r = (v — u)/2 we get equation (f2Ul) . 
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